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Произведены первопринципные расчеты характеристик взаимодействия водорода с алюминием. Исследовано влияние пере- 
стройки атомной структуры на величину энергии и объема растворения водорода в металле. В рамках модели химической свя- 
зи в металлах изучено влияние водорода на изменение электронной плотности в алюминии и проанализирован механизм воз- 
никновения избыточного объема, вносимого водородом в металл. 


1. Введение 

В связи с важностью описания механизмов изме- 
нения и прогнозирования свойств материалов в усло- 
виях взаимодействия с водородом [1-3] направлены 
усилия многих исследователей на изучение характе- 
ристик взаимодействия водорода с металлами. 

Энергия растворения водорода является одной 
из основных характеристик такого взаимодей- 
ствия, определяющая концентрацию примеси во- 
дорода при заданных внешних условиях и тем са- 
мым влияющая на все физические и механические 
свойства системы металл-водород. Избыточный 
объем, вносимый водородом в материал, предста- 


вляет другую величину, с которой связаны многие 
практические важные явления. Избыточный объем 
приводит к дальнодействующим взаимодействиям 
между атомами водорода и направленной мигра- 
ции водорода в поле упругих напряжений в разбав- 
ленных металл-водородных системах. Эти взаимо- 
действия ответственны за скопление водорода в 
областях с растягивающими напряжениями и вы- 
деление гидридов в полях напряжений вблизи вер- 
шин трещин, снятию дальнодействующих полей 
дислокаций и увеличения их подвижности. Рас- 
смотрение отмеченных и других явлений, связан- 
ных с энергией и объемом растворения водорода в 
материалах можно найти в книге [1] и обзорной ра- 
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боте [2]. В данной работе мы ограничиваемся изу- 
чением энергии и объема растворения одиночного 
атома водорода в ГЦК решетке алюминия на осно- 
ве расчетов характеристик системы водород-алю- 
миний из первых принципов. 

Первопринципные методы играют большую 
роль в получении детальной информации об осо- 
бенностях взаимодействия водорода с металлами в 
дополнение к результатам экспериментальных ис- 
следований. При этом большое число теоретиче- 
ских работ посвящено установлению особенностей 
взаимодействия одиночного атома водорода с ме- 
таллической матрицей как основы для детального 
понимания многих физических эффектов, имею- 
щих место в разбавленных металл-водородных си- 
стемах [2]. Алюминий занимает особое место в эт- 
их исследованиях в качестве модельного ГЦК ме- 
талла, достаточно простого для проведения расче- 
тов с большим количеством атомов в расчетной 
ячейке и в тоже время имеющий широкие техноло- 
гические приложения. 

Цель данной работы: дать физическую интер- 
претацию возникновения избыточного объема соз- 
даваемого водородом, до сих пор не рассматривав- 
шейся с привлечением первопринципных методов 
исследований, и рассмотреть зависимость между 
объемом и энергией растворения водорода в алю- 
минии. 

2. Метод расчета 

Расчеты рассмотренных в работе характеристик 
взаимодействия водорода с ГЦК решеткой алюми- 
ния проводились «из первых принципов» в модели 
сверхъячейки. На рис. 1 показаны положения те- 
тра- и октапор в ГЦК решетке, в центры которых 
помещались атомы водорода. Цифрами обозначе- 
ны атомы алюминия, через которые проходит пло- 
скость (111). Расчетная ячейка имела размеры 
2ох2ох2й в направлениях [100], [010] и [001], где 
а - постоянная решетки, и включала 32 атома алю- 
миния с одним атомом водорода в тетра- либо ок- 
тапоре. Выбранный размер ячейки достаточен для 
вычисления энергии растворения водорода в ГЦК 
решетке алюминия с точностью 0,01 эВ в методе 
сверхъячейки [4]. 

Первопринципные вычисления выполнялись, в 
рамках теории функционала плотности (ТФП) [5] с 
применением метода псевдопотенциала [6] и паке- 
та программ АВИЧІТ [7]. Обменно-корреляцион- 
ный функционал описывался в приближении ло- 
кальной плотности [8] с параметризацией предло- 
женной в работе [9]. Взаимодействие электронов с 
атомным ядром описывается псевдопотенциалом, 
имеющим форму Клеймана-Биландера [10]. Эта 
форма включает 5- и р-нелокальность и основана 
на первопринципном полулокальном псевдопо- 
тенциале Хаманна [11]. Полная энергия рассчиты- 
валась интегрированием по зоне Бриллюэна с ис- 
пользованием сетки точек 12x12x12, выбранной по 
схеме [12]. Заполнение электронных состояний 


Кона-Шема задавалось согласно схеме Ферми-Ди- 
рака с энергетическим параметром размытия по- 
верхности Ферми 0,025 эВ. Псевдоволновые функ- 
ции электронов задавились разложением по пло- 
ским волнам вплоть до энергии обрезания 700 эВ, 
что обеспечивало точность расчета энергии раство- 
рения водорода не хуже 0,01 эВ. Процедура самос- 
огласованного расчета электронной плотности 
проводилась до достижения точности в расчете 
полной энергии расчетной ячейки 2,7-10 5 эВ. Ре- 
лаксация атомной структуры к минимуму конфи- 
гурационной энергии принималась достигнутой, 
если силы, действующие на каждый атом, не пре- 
вышали 2,5- 10 3 эВ/А. При этом абсолютное значе- 
ние давления в системе не превышало ІО 3 Па, что 
определяло равновесный объем сверхъячейки. В 
результате расчетов получен равновесный пара- 
метр решетки для ГЦК А1 <7=3,96 А и объемный мо- 
дуль всестороннего сжатия 5=83,964 ГПа. Эти дан- 
ные находятся в хорошем согласии с эксперимен- 
тальными значениями а=4,05 А, 5=76,9 ГПа из ра- 
боты [14] и с результатами теоретических расчетов 
[13] (о=3,99 А и 5=83,2 ГПа). 
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Рис. 1 . Расположение междоузлий в ГЦК решетке АІ: ▲ - ок- 
таэдрическое и я ~ тетраэдрическое междоузлия, 
• - атомы АІ 

3. Результаты расчетов энергии 

и объема растворения водорода в алюминии 

Энергия растворения примеси водорода в ма- 
трице алюминия рассчитывалась по формуле 

А5 = 5(А1 32 )+І5(Н 2 )-5(А1 32 Н). (1) 

Здесь ЦАЦН) - полная энергия релаксирован- 
ной сверхъячейки с примесью водорода в тетра- ли- 
бо октапоре, ЦАЦ) - полная энергия сверхъячейки 
без примеси водорода и ЦН 2 ) - энергия двухатом- 
ной молекулы водорода, равная -30,88 эВ согласно 
нашим расчетам по изложенной выше методике. 
Формула ( 1 ) задает выигрыш в энергии при удале- 
нии атома водорода из молекулы и помещении его 
в междоузлие ГЦК решетки алюминия. Рассчитан- 
ные нами по этой формуле энергия растворения во- 
дорода в октаэдрическом междоузлии 
ДЦО)=— 0,53 эВ и в тетраэдрическом междоузлии 
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ДДТ)=— 0,46 эВ. Отрицательный знак энергии ра- 
створения указывает на термодинамическую неста- 
бильность примеси водорода в ГЦК решетке алю- 
миния в согласии с результатами недавних расчетов 
из первых принципов [4] и экспериментальными 
данными [15, 16]. При этом из наших расчетов сле- 
дует, что атому водорода энергетически выгодно за- 
нимать тетрапору. Разница энергий водорода в ок- 
таэдрическом и тетраэдрическом междоузлиях рав- 
на 0,07 эВ, что хорошо согласуется с соответствую- 
щей величиной 0, 1 3 эВ, полученной авторами рабо- 
ты [4]. Небольшое количественное расхождение 
можно объяснить использованием в работе [4] при- 
ближения сверхмягкого псевдопотенциала. 

В результате растворения атома водорода в кри- 
сталле алюминия появляется избыточный объем, 
который в методе сверхъячейки определяется раз- 
ницей объемов сверхъячеек с примесью водорода и 
в отсутствии примеси. При этом расчетная величи- 
на избыточного объема не зависит от выбора свер- 
хъячейки, если она имеет достаточно большой раз- 
мер. В работе [4] было установлено, что расчетная 
величина энергии растворения водорода практиче- 
ски не меняется (в пределах 0,01 эВ) для сверхъяче- 
ек с 32, 64 и 108 атомами алюминия. Это указывает 
на адекватность описания релаксацинных смеще- 
ний атомов уже в сверхъячейке из 32 атомов алюми- 
ния. Поэтому мы проводили расчеты для сверхъя- 
чейки этого размера и, можно утверждать, что рас- 
считанные нами величины избыточных объемов 
4,81 и 3,22 А 3 для водорода в тетраэдрическом и ок- 
таэдрическом междоузлиях, соответственно, каче- 
ственно правильно отражают изменения объема, 
вносимые водородом в алюминий. Для сравнения 
приведем величину избыточного объема 3,97 А, ко- 
торую вносит водород в октаэдрическом междоуз- 
лии алюминия, полученную авторами работы [17] 
из расчетов методом погруженного атома. Результа- 
ты других расчетов избыточного объема водорода в 
металлах из первых принципов нам не известны. 
Экспериментальные данные по избыточному объе- 
му водорода в алюминии также отсутствуют из-за 
его слабой растворимости в этом металле, что не 
позволяет провести измерения в рамках суще- 
ствующих экспериментальных методик [17]. Ре- 
зультаты наших расчетов качественно согласуются 
с экспериментально определенными избыточными 
объемами водорода в других металлах: 2,48 А (N1), 
2,78 А (Ті), 3,11 А (2г) и 2,9 А (Ре) (данные взяты 
из работы [17]). Большее значение рассчитанной 
нами величины избыточного объема коррелирует с 
меньшей растворимостью водорода в алюминии в 
сравнении с перечисленными металлами. 

Происхождение вносимого водородом избы- 
точного объема с точки зрения химической связи в 
металлах можно понять из исследования измене- 
ний зарядовой электронной плотности, вызванных 
растворением водорода в металлах. 

На рис. 2 изолиниями показана разность вели- 
чин электронных плотностей в плоскости (111) в 


сверхъячейках с атомом водорода и без примеси во- 
дорода при идентичных положениях атомов алю- 
миния в узлах ГЦК решетки. Цифрами обозначены 
три ближайшие к атому Н (помещенного в тетрапо- 
ру) атома А1, через которые проходит плоскость 
(111), показанная ранее на рис. 1. Таким образом, из 
рис. 2 видно, что вносимая электронная плотность 
резко возрастает с приближением к протону (кон- 
турные линии в этой области располагаются на- 
столько близко, что сливаются в черное пятно; что- 
бы не загромождать рисунок этой деталью мы соч- 
ли за необходимость представить его белым пят- 
ном), что связано с экранировкой его положитель- 
ного заряда. В то же время наблюдается увеличение 
электронной плотности также и в соседних меж- 
доузельных областях, не содержащих атом водоро- 
да. Аналогичное распределение вносимой водоро- 
дом электронной плотности наблюдалось нами и в 
других плоскостях возле атома водорода, в том чи- 
сле и для атома водорода в октаэдрическом меж- 
доузлии. Результаты соответствующих расчетов из- 
менений электронной плотности здесь не приво- 
дятся из-за их качественного подобия результатам, 
представленным на рис. 2. Таким образом, харак- 
терной особенностью изменения электронной 
плотности, связанной с внесением атома водорода, 
является ее возрастание в междоузельных областях. 



Рис. 2. Разность зарядовых плотностей валентных электро- 
нов. Контуры плотности даны в единицах е/П с мас- 
штабом 0,001, где е - заряд электрона, а О. ~ объем 
элементарной ячейки, а • , ▲ - положение атомов АІ 
и Н соответственно 

Поведение электронной плотности в меж- 
доузельных областях составляет важный элемент 
установления равновесного объема металлов в рам- 
ках модели металлической связи, развитой на осно- 
ве многочисленных расчетов электронной структу- 
ры [18]. Согласно этой модели выигрыш в зонной 
энергии электронов обуславливает стремление ато- 
мов металла к уменьшению расстояний между ни- 
ми, что приводит к увеличению электронной плот- 
ности в междоузельных областях и возрастанию да- 
вления электронного газа, зажатого в междоузлиях. 
В результате противодействия этих тенденций уста- 
навливается равновесный объем металла при неко- 
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торой величине междоузельной электронной плот- 
ности [18]. Развитая модель успешно объяснила на- 
блюдаемые закономерности в различиях величин 
модуля всестороннего сжатия металлов, которые 
хорошо согласуются со значениями модуля всесто- 
роннего сжатия однородного электронного газа с 
плотностью, равной средней электронной плотно- 
сти в междоузлиях металла [18]. 

Увеличение электронной плотности в междоуз- 
лиях при внесении примеси водорода нарушает су- 
ществующий баланс между сжимающими и растяги- 
вающими тенденциями в металле. Для установления 
нового равновесного объема решетка расширяется в 
области воздействия атома водорода под действием 
увеличенного давления со стороны электронного га- 
за в междоузлиях с более высокой средней электрон- 
ной плотностью. Это объясняет появление избыточ- 
ного объема, вносимого водородом, в рамках рассмо- 
тренной выше модели химической связи в металлах. 
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Рис. 3. Смещение атомов алюминия вдоль кристаллографи- 
ческих направлений, где А = ((К- Ко)/ Ко)- 100%, а 

Ко, К - положение атомов АІ относительно атома Н до 
и после релаксации атомной структуры соответственно 

Интересно отметить, что согласно нашим расче- 
там абсолютная величина энергии растворения в 
тетраэдрическом междоузлии (-0,46 эВ) меньше, 
чем в октаэдрическом междоузлии (-0,53 эВ), в то 
же время величины объемов растворения в тетрапо- 
ре (4,81 А 3 ) больше чем в октапоре (3,22 А). Отме- 
ченная особенность требует объяснения на атом- 
ном уровне, поскольку в континуальной среде, не- 
сомненно, существует однозначная связь - чем 
больше объем образования точечного дефекта, тем 
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больше его энергия образования. Рассмотрим осо- 
бенности релаксационных смещений атомов во- 
круг атома водорода, расположенного в тетраэдри- 
ческом и октаэдрическом междоузлиях. На рис. 3 
представлены величины смещений атомов алюми- 
ния по кристаллографическим направлениям от 
атома водорода. Как видно из рисунка, смещения 
от атома водорода в тетраэдрическом междоузлии 
больше по величине и распространяются на боль- 
шие расстояния, что объясняет происхождение 
большего избыточного объема, вносимого атомом 
водорода при занятии тетраэдрического междоуз- 
лия. С другой стороны, большие смещения атомов 
указывают и на более значительное уменьшение 
энергии системы в результате релаксации атомной 
структуры. Представленный анализ смещений ато- 
мов согласуется с результатами расчетов энергий 
релаксации атомной структуры после помещения 
атома водорода в междоузлия (разницы энергии 
сверхъячейки при расположении атомов алюминия 
в узлах идеальной решетки с атомом водорода в 
междоузлии и энергии той же сверхъячейки с рав- 
новесными положениями атомов) [18]. Без учета 
энергии релаксации разница энергий водорода в те- 
траэдрическом и октаэдрическом междоузлиях по- 
ложительна, а в результате атомной релаксации эта 
разница меняет знак. Таким образом, однозначная 
связь между энергией и объемом растворения водо- 
рода в металле отсутствует, и эти параметры опреде- 
ляются особенностями в перестройке атомной 
структуры, вызванной вносимым водородом. 

4. Заключение 

Расчеты из первых принципов энергии раство- 
рения показывают, что водород занимает тетра- 
эдрическое междоузлие в ГЦК решетке алюминия 
и является термодинамически нестабильной при- 
месью в кристалле. Обнаружена корреляция в ве- 
личинах энергии и объема растворения Н в АІ. 
Установлено, что такие особенности обусловлены 
значительно большими смещениями атомов ме- 
талла при занятии водородом тетраэдрического 
междоузлия. Из анализа изменений самосогласо- 
ванной электронной плотности, вносимых атомом 
водорода в алюминий, дана физическая интерпре- 
тация появления избыточного объема в металле в 
рамках модели химической связи в металлах. 
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